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Abstract. We study the dynamics of solitons as solutions to the perturbed KdV 
(pKdV) equation dtu = —dxid^u + — bu), where b{x,t) = bo{hx,ht), ft, <C 1 
is a slowly varying, but not small, potential. We obtain an explicit description of 
the trajectory of the soliton parameters of scale and position on the dynamically 
relevant time scale 6h~^ log , together with an estimate on the error of size 
h^/^. In addition to the Lyapunov analysis commonly applied to these problems, 
we use a local virial estimate due to Martel-Merle [12]. The results are supported 
by numerics. The proof does not rely on the inverse scattering machinery and is 
expected to carry through for the subcritical gKdV-p equation, 1 < p < 5. The 
case of p = 3, the modified Kortcwcg-dc Vries (mKdV) equation, is structurally 
simpler and more precise results can be obtained by the method of Holmer-Zworski 



1. Introduction 
The Korteweg-de Vries (KdV) equation 

(1.1) dtu = d^{-dlu - 3u'^) 

is globally well-posed in if ^' for A; > 1 (see Kenig-Ponce-Vega [13]). It possesses soliton 
solutions u{t, x) = r]{x, a+Ac^t, c), where r]{x, a, c) = c^9{c{x—a)) and 9{y) = 2 sech^ y 
(so that 6" + 36^ = 46). Benjamin [1], Bona [2], and Bona- Souganidis- Strauss [3] 
showed that these solitons are orbitally stable under perturbations of the initial data. 
We consider here the behavior of these solitons under structural perturbations, i.e. 
Hamiltonian perturbations of the equation (II. ip itself. Dejak-Sigal [4], motivated 
by a model of shallow water wave propagation over a slowly-varying bottom, have 
considered the perturbed KdV (pKdV) 

(1.2) dtU = d^{-dlu - 3u^ + bu) 

where b{x,t) = h^~^^bo{hx, ht) and h <^ 1. They proved that the effects of this 

potential are small on the dynamically relevant time frame. We consider instead 
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b{x,t) = bQ{hx,ht), a slowly- varying but not small potential jl| which allows for con- 
siderably richer dynamics. 

To state our main theorem, we need the following definition: 

Definition 1 (Asymptotic time-scale). Given b^ E C^(M?), G M, Cq > 0, and 

6 > 0, let A[t), C{t) solve the system of ODEs 

ri = 4c2-6o(A-) 

(1-3) { . 1 

with initial data A(0) = Aq and C(0) = Cq. Let T* be the maximal time such that on 
[0,T^,), we have 5 <C{t) < 5^^ . (T^ could be +oo.) 

Let (m, v) = f uv. 

Theorem 2. Given bo G C^{M.'^), Aq G M, Co > 0, and 6 > 0, let T, be the 
time defined in Def. [H Let oq = h~^Ao and cq = Co- Then for < t < T =^ 
min(T,,, 51og /i~^) , there exist trajectories a{t) and c{t), and positive constants 
e = e{6) and C = C{6,bo), such that the following holds. Taking u{t) the solution 
of ( II. 2p with potential b{x,t) = bo{hx,ht) and initial data ri{-,ao,Co) , let v{x,t) =^ 
u{x,t) — ri{x,a(t),c(t)). Then 

(1-5) \\e-'^''-MLf,^,,Hi<h'^'e^'\ 
and 

(1.6) (w,?7(-,a,c)) = 0, {v,{x- a)7]{-,a,c)) = 0. 
Moreover, 

(1.7) \a{t) - h~'A{ht) I < e^^* , |c(t) - C(/it) | < he^'''' . 

Up to time 0{h~^), a(t) is of size 0{h~^) and c{t) is of size 0(1), and (11.71) gives 
leading-order in h estimates for a{t) and c{t) - that is, despite the differences in 
magnitudes, the estimates for a{t) and c{t) provided by (11. 7p are equally strong. The 
strength of the local estimate (11.51) . in comparison to the global estimate (ll.4p on the 
error v, is that it involves integration in time over a (long) interval of length 0{h~^). 
The estimate (11.50 is on par, although slightly weaker than, the pointwise-in-time 
estimate ||e~''''^'~'*'f 2,2 < /le'"'**. The two estimates ( 11. 4p . ( 11. 5p are consistent (but 

[0,T] a: 

not equivalent to) v being of amplitude h but effectively supported over an interval of 



""^Dejak-Sigal [3| state a more general result that appears to allow for potentials that are not 
small. However, the smallncss in their result is required to reach the dynamically relevant time 
frame ^ . See the comments below in fJHH 
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size 0{h~^), which is suggested by numerical simulations. The trajectory estimates 
( 11. 7p state that we can predict the center of the soliton to within accuracy 0(1) and 
the amplitude to within accuracy 0{h). (This discussion does not include the h^^ 
loss that occurs when passing to the natural Ehrenfest time scale 5h~^ \ogh~^ .) 
To define the Hamiltonian structure associated with (11. 2p . let J = dx with 

j-V(x) = 9-V(x) = - / -/ ]fiy)dy. 



We regard the function space = H^{M.) as a symplectic manifold with symplectic 
form u{u, v) = {u, J~^v) densely defined on the tangent space TN ~ H^. Then (11. 2p 
is the Hamilton fiow dtu = JH'{u) associated with the Hamiltonian 

(1.8) ^ ^ \ ^ {ul-2u' + hu' 



X 



Let M G N = denote the two-dimensional submanifold of solitons 

M = {r/(-,a,c) I a e M, c> 0} . 

By direct computation, we compute the restricted symplectic form uj\^j = Sc^da A dc 
(thus M is a symplectic submanifold of TV) and restricted Hamiltonian = — yc^ + 
^B{a, c, t), where 

B{a, c, t) =^ J b{x, t)ri{x, a, c)^ dx . 

The heuristic adopted in [8l|9], essentially equivalent (see [10]) to the "effective La- 
grangian" or "collective coordinate method" commonly applied in the physics litera- 
ture, is the following: the equations of motion for a, c are approximately the Hamilton 
flow of if I ^ with respect to These equations are 

16 

: = —C-''daB 

16 

By Taylor expansion, these equations are approximately 

= - b{a) + 0{h^) 

= -cb'{a) + 0{h^) 
3 

Note that the equations (11.31) are the rescaled versions of these equations with the 
0{h?) and 0{h^) error terms dropped. 

The flrst of the orthogonality conditions in (11.61) can be rewritten as a;(w, daTj) = 
and thus interpreted as symplectic orthogonality with respect to the a-direction on 
M. The other symplectic orthogonality condition = uj{v, dj]) = {v, d'^dcf]) is 
not deflned for general functions v since d~^dcr] = {T{y) + y9{y))\ where 

\y — c[x a) 

T{y) = 2 tanhy. Thus, we drop this condition, although it must be replaced with some 
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other condition that projects sufficiently far away from the kernel (spanl^^.//}) of the 
Hessian of the Lyapunov functional. We select (f , {x — ajrj) = (i.e., the second 
equation in (11. 6p ) since it is a hypothesis in the Martel-Merle local virial identity 
(Lemma I6.ip . 

1.1. Numerics. For the numerics, we restrict to time- independent potentials b{x) = 
bo{hx) and use the rescaled frame X = hx, S = hH, V{X,S) = h~^u{h~^X, h~^S), 
and B{X) = h-%{h-^X) = h-%o{X). Then V solves the equation 



with initial data Vq{X) = ri{X, Aq, CqH Note that to examine the solution u(x, t) 
on the time interval < t < Kh~^, we should examine V{X, S) on the time interval 



As an example, we put bo{x) = 8 sinx and take Aq = 2.5, Cq = 1 and K = 1. Then 
the width of the soliton is approximately the same width as the potential (when 
h = 1), but note that the size of the potential is not small. The results of numerical 
simulations for h = 0.3,0.2,0.1 are depicted in the Fig. [TJ There, plots are given 
depicting the rescaled solution v{X, S) for each of these values of h. In Fig. [2], we 
draw a comparison to the ODEs (II. 3p . In each of the numerical simulations, we record 
the center of the soliton as Ah{S) and the soliton scale as 



That is, we ffi the solution V{X,S) to ri{X, Ah{S),Ch{S)). Let T = /it so that 
S = h?T . To convert into the {X^T) frame of reference, we plot T versus Ah{T) = 
Ah{h^T) in the top plot of Fig. [2] together with A{T) solving ([L3]). In the bottom 
frame, we plot T versus Ch(T) = hCh{h?T) together with C{T) solving (11.31) . We 
opted to only plot h = 0.2 since the curves for h = 0.3, 0.2, 0.1 were all rather close, 
producing a crowded figure. Theorem [2] predicts 0{h) convergence in both frames of 



The numerical solution to the equation (ll.2p was produced using a MATLAB code 
based on the Fourier spectral/ETDRK4 scheme as presented in Kassam-Trefethen 
|llj . The ODEs (II. 3p were solved numerically using 0DE45 in MATLAB. 

1.2. Relation to earlier and concurrent work. Theorem 2 in Dejak-Sigal [1] 
states (roughly) that for potential b{x,t) = eb{hx, ht), the error < e^l'^h}!'^ can 

be achieved on the time-scale t < (/i-|-e^/^/i^/^)~^, and the equations of motion satisfy 



To reach the nontrivial dynamical time frame, one thus needs to take e = /i in their 
result. With this selection for e, the 0(/i^) errors in the ODEs can be removed as 



dsV = dxi-d^V - 3V^ + BV) , 



< 5 < Kh^. 




Fig. H 



d = - b{a) + 0{eh) 
c = 0{eh) 
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FlGURE 1. The rescaled evolution V{X,S) (see text) for B{X) = 
h~%Q{X) = 8h~^smX, Aq = 2.5, Cq = 1, on the time interval 
< S < h^. The three frames are, respectively, h = 0.3, h = 0.2, 
and h = 0.1. 
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Figure 2. For the simulations in Fig. [H the position was recorded 
as Ah{S) and the scale was recorded as Ch{S); that is, the solution 
v{X,S) was fitted to r]{X, Ah{S),Ch{S)). The top plot is T versus 
Ah(T) = Ah{h?T) for h = 0.2 (in blue) compared to the value of A(T) 
obtained by solving the ODE system (in green). The bottom plot is T 
versus Ch(T) = hCh{h?T) for h = 0.2 (in blue), compared to the value 
of C(T) obtained from the ODE system (in green). 
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in our result with the effect of at least preserving the error estimate for w in at 
the /i^/^, rather than h level. But then the conclusion of their analysis is that the 
(small and slowly varying) potential has no significant effect on the dynamics. We 
emphasize that in our case, we allow for e = 0(1) and thus can see dramatic effects 
on the motion of the soliton. 

The paper Dejak-Sigal [1] is modeled upon earlier work by Frohlich-Gustafson- 
Jonsson-Sigal [5] for the NLS equation, which controlled the error via the Lyapunov 
functional employed in the orbital stability theory of Weinstein [18] . In |9] , we im- 
proved |5] by using the symplectic restriction interpretation as a guide in the analysis 
and introducing a correction term to the Lyapunov estimate. A correction term is not 
as easily applied to the study of (11. 2p since the leading order inhomogeneity in the 
equation for v generates a "nonlocal" solution. To properly address the nonlocality 
of V, we use both the global estimate (ll.4p as in [SI HI [HI E] , but also introduce 
the new local estimate (II. 5p . which is proved using the local virial identity of Martel- 
Merle [15]. We remark that our method does not use the integrable structure of the 
KdV equation, and we expect that our result will carry over to the perturbed 
subcritical gKdV-p equation 

dtu = —dx{d1u + vF — hu) 

In the case p = 3, i.e. the second symplectic orthogonality condition (f , d'^dcfj) = 
(where now 9(y) = \/2sech.y and r](x, a, c) = cr]{c{x — a))) is well-defined for general 
functions v. In this case, we are able to achieve stronger results by following the 
method of [9], and even treat double solitons - see [7]. 

The concurrent work by Munoz [T7] considers the equation (specializing to the case 
m = 2 in his paper to facilitate comparison) 

(1.9) dtv = d^i-d^v + AXv - 3av^) , 

where a{x) = aoi^hx), with ao(^) increasing monotonically from a{—oo) = 1 to 
a(oo) = 2, and < A < Aq == | is constant, effectively corresponding to a moving 
frame of reference. The equation (II. 9p is similar to our (II. 2p but not directly related 
to it through any known transformation. His main theorem gives the existence of 
a solution v{x,t) which asymptotically matches the soliton ri{x,a{t), 1) as t — )■ — oo 
and matches the soliton ^r]{x,a{t), Coo) as t — j- +oo with error at most h^^'^ in if^. 
Here, Coo is precisely given in terms of the solution to an algebraic equation (see (4.17) 
in his paper). He presents this problem as more of an obstacle scattering problem 
with a careful analysis of "incoming" and "outgoing" waves and thus his priorities 
are different from ours. 

However, information from the "interaction phase" of his analysis can be extracted 
from the main body of his paper and compared with the results of our paper. In the 
course of his analysis, he obtains effective dynamics (here Aq = |) for an approximate 
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solution 

{a = (? — \ 
5 V -^0/ 

He then shows that the approximate solution is comparable to a true solution in 
with accuracy 0{h}/'^) (same as in our result) but only at the expense of a spatial 
shift for which he has the comparatively weak control of size 0{h~^). In our analysis, 
we are able to achieve control of size 0(1) on the positional parameter a{t). At 
the technical level, we are gaining an advantage by using the local virial estimate in 
the interaction phase analysis while Muiioz carries out a more direct energy estimate. 
Mufioz does apply the local virial estimate in his "post-interaction" analysis to achieve 
a convergence statement as t — )■ +oo with a remarkably precise scale estimate. 

1.3. Notation. It is convenient to work in both direct (e.g. ?7(x, a,c)) and "pulled- 
back" coordinates (e.g 9{y)). Our convention is that successive letters are used to 
define functions related in this way. Specifically, 

• 6{ii) = 2 sechy and ri{x, a, c) = c^9{c{x — a)). 

• '''{y) = 2 tanhy and a{x, a, c) = c^r(c(x — a)). 

• v{x, t) = 2c^w{c{x — a),t) 

• £ = 4 - 92 - 6^ and /C = 4c2 - 92 - 6r/(-, a, c). 

1.4. Outline of the paper. In ^ we deduce some needed spectral properties of the 
operator /C which are required to give the lower bound in the Lyapunov functional 
method (Cor. 12. 4p . In ^ we give the standard argument, via the implicit function 
theorem, that the parameters a and c can be adjusted so as to arrange that v satisfies 
the orthogonality conditions ( II. 6p (Lemma 13. ip . In ^ we decompose the forcing term 
in the linearized equation into symplectically orthogonal and symplectically parallel 
components. In ^ the orthogonality conditions are applied to obtain the equations 
for the parameters (Lemma 15. ip . These equations include error terms expressed in 
terms of the local-in-space norm ||e~'^'^~"'f In ^ an estimate on ||e~'^'^~"'w||/,2 j^^i 
is obtained by the Martel-Merle local virial identity (Lemma 16. 3p . In ^ the estimates 
on llt'lli^j/i are obtained by the Lyapunov energy method (Lemma 17. ip . The three 
key estimates (Lemmas 15. H 16. 3[ 17. ip are combined to give the proof of Theorem [2] in 

1.5. Acknowledgements. Galina Perelman shared with me a set of notes illustrat- 
ing how to apply the Martel-Merle local virial identity to this problem. The present 
paper is essentially an elaboration of this note, and hence I am very much indebted 
to her generous assistance. I thank also Maciej Zworski for initially proposing the 
problem, providing the numerical codes, and for helpful discussions. 

I am partially supported by a Sloan fellowship and NSF grant DMS-0901582. 
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2. Spectral properties of the linearized operator 

Recall that C = A — dy — 69. Since 6{y) = 2 sech^ y, we see that we must consider 
the Schrodinger operator with Poschl- Teller potential 

A = -dy - u{u + 1) sech^ y 

with z/ = 3. The spectral resolution of operators of the type A is deduced via hyper- 
geometric functions in the appendix of Guillope-Zworski [6]. From this analysis, we 
obtain 

Lemma 2.1 (spectrum of £). The spectrum of C is {—5,0,3} U [4, +oo). The 
normalized eigenfunctions corresponding to the first two eigenvalues are 



Ai = -5 fi{y) = sech^ y 



Ao = fo{y) = -^sech^ ytanhy = ^0'{y) 

z o 

Denote by Ej the corresponding eigenspaces and Pej the corresponding projections 
(that is, the orthogonal projections and not the symplectic orthogonal projections). 

Lemma 2.2. Suppose that {w,6) = and {w,y9) = 0. Then 

(2.1) 2\\w\\l2 < {Cw,w) 

Proof. Since C preserves parity, it suffices to separately prove: 
Claim 1. If w is even, ||l2 = 1, and {w,6) = 0, then {Cw,w) > 2. 
Claim 2. If w is odd, ||u^||l2 = 1, and {w,y6) = 0, then {Cw,w) > 2. 
We begin with the proof of Claim 1. Since w is even, {w, fo) =0. Resolve w as 

w = afi + g, ge{Ei + Eo)^, + \\g\\l, = 1 . 

Resolve also 



e = f3fi + h, he{Ei + Eo)^, P'+\\h\\h 
We compute that 



2 16 



L2 



3 



3a/1 ^tt 

(2.2) /3 = {6 Ji) = ^-^^^ 2.28138, 



16 



from which it follows that 



(2.3) ll^lli. = f - j - 0.128659. 

We then have 

= {w,9)=a/3+{g,h), 



10 JUSTIN HOLMER 

which using f l2.2p . f l2.3p . and ||5'||l2 < 1, imphes 

|a| < ^IklUHI^IlL^ < 0.157226. 

By the spectral theorem, 

w) > 3\\g\\l2 - = 3(1 - a^) - 5a^ = 3 ~ 8a^ > 2 . 
Next, we prove Claim 2. Since w is odd, {w, fi) = 0. Resolve w as 
w = afo + g, ge{Ei + Eo)-^, a"^ + \\g\\l2 = 1 ■ 

Resolve also 

ye = Pfo + h, he{E,+E,)^, p^ + \\h\\i, = \\ye\\l, = ^{Ti 
We compute that 

(2.4) /5 = (t/0,/o) = 'y/|^ 1.29099, 
from which it follows that 

(2.5) \\h\\l2 = -(vr^ -Q)-f3^^ 0.0531575 . 

9 

We then have 

0= {w,ye) = al3+{g,h) , 
which, using (12. 4p . (12. 5p . and \\g\\L'^ < 1 implies 

|«| < ^II^IUHI^iU^ < 0.17859. 

By the spectral theorem, 

w) > 3\\g\\l2 = 3 - 3a2 > 2 . 



Corollary 2.3. Suppose that 

(2.6) {w,9) = and {w,y9) = 

Then 

11 1 

Proof. By integration by parts. 



{Cw,w) = 4:\\w\\l2 + \\dM\h - G J Ow^ 



from which we obtain 



\dM\h < {Cw,w) + 8\\wf, 



lL2 

Adding to this estimate |x the estimate (12. ip . we obtain the claim. 
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Of course the above properties of £ can be converted to properties of /C, where 

)C = Ac^-dl-Q7]{-,a,c), 
by scahng and translation. In particular, we have 
Corollary 2.4. Suppose that 

(2.7) {v,ri{-,a,c)) = and {v, {x — a)ri{-,a,c)) = . 

Then 

\\v\\hi < {}Cv,v) . 
where the implicit constant depends on c. 

3. Orthogonality conditions 

We next show by a standard argument that the parameters (a, c) can be tweaked 
to achieve the orthogonality conditions (11. 6p . 

Lemma 3.1. If S < c < d^^, there exist constants e > 0, C > such that the 
following holds. If u = r]{-,d,c) +v with \\v\\h^ < e, then there exist unique a, c such 
that 

\a — d\ < C\\v\\Hi , \c — c\ < C\\v\\h'^ 
and f =^ M — T]{-,a, c) satisfies 

(f,?7)=0 and {v, {x — a)ri) = . 
Proof Define a map $ : if ^ x M x M+ by 

$(« ac)=\ (^-^(■'«'C),r^) 

[{u- ri{-,a,c),{x - a)T]) 

The derivative of $ with respect to (a, c) at the point (?7(-, a, c), a, c) is 

{daV, v) {dcV, V) 

{daT], (x - a)r]) {dcT], (x - a)r]) 

which is nondegenerate. By the implicit function theorem, the equation $(u, a, c) = 
can be solved for (a, c) in terms of m in a neighborhood of r]{-, d,c). □ 

4. Decomposition of the flow 

Since we will model u = ri{-, a,c) + v and u solves (11.21) . we compute that v solves 
dtv = —dx{d1v + Qr]v — hv + 3f^) + Fq 
(4.1) = dXv - Ac^d.,v + dx{hv) - W^v^ + Fo 

where 

Fo = -(d - 4:C^)da'n - cdj] + 9^ (677) . 



{Da,c^){ri{-,a,c),a,c) 





' 8c2" 




.V _ 
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Decompose Fq = F\\ + F±, where F\\ is symplectically parallel to M and F± is sym- 
plectically orthogonal to M. Explicitly, we have 

By Taylor expansion we obtain Fj_ = (-Fl)o + 0{h?), where 

{F^)o = \c'y{a) (%) + 2y^'(y))|^^^(_^. 

By definition of F^, we have {F^,d~^dar]) = and {F^,d~^dcrj) = , which must 
then hold at every order in h; in particular, they hold for (-Fl)o- Note that by parity 
(Fj_)o in addition satisfies ((-Fj_)o, {x — a)ri) = 0, although this is not expected to hold 
for F± at all orders. 
It follows that 

(4.2) ||e^l^-"lFo|ki <\a- ^c' - b{a)\ + \c - ic6'(a)| + h . 

5. Equations for the parameters 

Lemma 5.1. Suppose that we are given G C^(]R^) and 6 > 0. (Implicit constants 
below depend only on and 6.) Suppose that \\v\\h]^ 1, v solves (14.11) and satisfies 
^;6^ , and 6 < c < 6-\ Then 

(5.1) \c-^cb'{a) \ < h\\e-'^^-''^v\\Hi + ||e"^l^"'*lt;||^i + /i^ 

and 

{dr^JCv, {x - a)ri) 



(5.2) 



d - 4c^ + 6(a) + 



{d^7], {x - a)r]) 



Proof. We first work with the orthogonality condition (f , ^daf]) = to obtain (15. ip . 
Applying dt to this orthogonality condition, we obtain 

= {dtv, ri{-, a, c)) + {v, dtr]{-, a, c)) . 
Substituting the equation for v and the relation dtt] = adaf] + cdcf], we obtain 
= {d^lCv, T]) - Ac^{d,v, v) + {d.M^v) - Hd.y, r]) ^ I + II + HI + IV 

+ r/) + (Fx, v) + a{v, daV) + c{v, d^v) ^ V + VI + VII + VIII 

We have 1 = and II = 0. Next, we calculate 

III = {d,{bv),rj) = -{bv,ri') = b{a){v,ri') + 0{h)\\e-^\--Mm 
= 0{h)\\e-^^^-MH^ 
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We easily obtain |IV| < ||e"'l^"''lw||^i. Next, 



V=(F|,,r/)=(-c+^9,i?)(f},r/,r/) 



from which it follows that 



y = -Sc\c - Icb'ia)) + 0{h^) . 
Next, we have VI = and Vll = 0. Finally, 

|V111| < \c-lcb'ia)\\\e^^\^-MHi + h\\e^'\^^MHi. 

Using that ||f Hi^i ^ 1, we obtain flS.ip . 

To establish fl5.2l) . we apply dt to {v, (x — a)r]) = to obtain 

= {dtv, {x - a)ri) + {v, dt[{x - a)rf\) 

Substituting the equation fl4.1l) for v and the relation dtf] = adat] + cdcTj, we obtain 

= {d^lCv, (x - a)r]) - Ac^{d^v, (x - a)r]) + {d^{hv), (x - a)r]) ^ 1 + II + III 
- 3(9^^;^ (x - a)ri) + {F\\ , (x - a)r/) + (F^, (x - a)r/) ^ IV + V + VI 

+ a{v, da[{x - 0)7]]) + c(f , (x - a)dcr]) ^ VII + VIII 

Note that we do not have 1 = 0. We would have I = if we were working with the 
orthogonality condition (w, d'^dcTj) = 0, but as explained previously, this condition 
cannot be imposed on v via the method of Lemma \37L\ and even if it could, would not 
give the coercivity in Corollary 12.41 We therefore keep Term I as is for now. Next, 
we note that 

II + III + VII = (-4c2 + b{a) + d){d^v, (x - a)r]) + 0{h)\\e-'\''-%\\Hi ■ 
Next, |IV| < ||e-^l^-"li;||^i. Also, 

V=(^-a + 4c' - j^d.B^ {daV, {x - a)r/) . 

It happens that {6 + 2y9',ye) = and hence VI = 0{h'). Finally, |VIII| < 
|c|||e~^'^~"'f 11/^1, to which we can append the estimate (15.11) . Collecting, we obtain 

□ 
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6. Local virial estimate 

Next, we begin to implement the Martel- Merle [15] virial identity. Let $ G C(M), 
<|)(x) = < on (0, +oo) such that = 1 on [0, 1] and $(x) = e"^ on 

[2, +oo), and e"^ < <l>{x) < Se"^ on (0, +oo). Let ^{x) = <!>{y) dy, and for A > 1 
(to be chosen later) set iplx) = A'^lx/A). 

The following is the (scaled-out to unity version of) Martel-Merle's virial estimate. 

Lemma 6.1 (Martel-Merle [IHl Lemma 1, Step 2 in Apx. B] and [HI Prop. 6]). 

There exists A sufficiently large and Aq > sufficiently small such that if w satisfies 
the orthogonality conditions 

{w,e) = and {w,ye) = 0, 

then we have the estimate 



An 



Step 2 in Apx. B of [15] is a localization argument that shows that it suffices to 
consider the case A = oo and ip{y) = y. Some integration by parts manipulations 
and the fact that {w, 6) = convert this case to the estimate 

(6.1) Wwfm < hlw, w) + j^{w, y9'){w, d') , 

where L = (| + 2y6' — 26) — dy. The positivity estimate ( 16. ip appears as Prop. 3 in 
[T5] and as Prop. 6 in [13], and is proved in [T^ . 
By scaling Lemma [6. H we obtain the following version adapted to /C. 

Corollary 6.2. There exists A sufficiently large and Aq > sufficiently small such 
that if V satisfies the orthogonality conditions (11. 6p . then (with tp = ■il){x — a)) 

Lemma 6.3 (application of local virial identity). Suppose that we are given bo G 
C^{M?) and 6 > 0. (Implicit constants below depend only on bo and 6.) Suppose that 
I — d + 4c^ — 6(a)| ^ 1, \\v\\ui <^ 1, V solves (14. ip and satisfies (II. 6p . and S < c < 6~^. 
Then with ip = ip {x — a) , we have 

(6.2) \\e~'^"~M\m^ < -f^idt j + ^^2h'' + H2h\\v\\l, , 

where e = e((5) > and j((5, 6o) > 0. Integrating over [0,T], we obtain with 

T < 

(6.3) \\e"^'"M\L^ m<\HL^ m+T^'^h. 
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Proof. Recalling that = iplx — a), 

dt J 'ipv^ = -a j i)'v'^ + 2 J ipvd^Kv - Sc^ J ipvd^v + 2 / tpvd^{bv) 

We reorganize the terms in the equation to 
-2 / / ^vFq 



A B 







I II in IV V 

Note that we have written this equation symbolically in the form 

(6.4) A + B = I + II + III + IV + V, 

and we now consider these terms separately. Integration by parts yields 

III = / ^'w^ 



IV = - / i)'bv'^ + j tpKv'^ 

ip'b{a)v^ - [ tp'{b{x) - b{a))v^ + / tpb^v^ 



Hence 



II + III + IV = (-d + - 6(a)) J ^'v^ + 0{h)\\v\\l2 , 



from which it follows that 

(6.5) |II + III + IV| < |d-4c2 + 6(a)|||e-^l"-"lt;||2, 
Integration by parts also yields 

V = 4 y ^'v^ , 

from which it follows that 

(6.6) \y\ < We^'^^'^MliMm. ■ 
Using that 

Fo = (d - 4c2 + b{a))dxri + 0(/i + |c|)e-2^l"-'^l , 

we obtain 

B = -2(d - 4c' + b{a)){ijv, dxT,) + 0{h + \c\)\\e-'\^-MLl 
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By (O, 

(6.7) B ^ 2i^^^(,., OiH + |c|)||e-l-^l.|U, . 

Placing estimates fl6.5p . (16. 6p . and (16. 7p into (16. 4p . we obtain, for some constant k > 0, 
the bound 

{d^f], [x - a)ri) 

< -9tJ iJv^ + K{\d-Ac' + b{a)\ + \\v\\Hi)\\e-'^''-M\k 

+ K{h+\c\)\\e-'^''-''^v\\L2+Kh\\v\\l2 

Using Corollary 16.21 and the assumptions |d — 4c^ + 6(a) | ^ 1, ^ 1, we obtain, 

for some constants Hi, H2 > 0, the bound 

(6.8) \\e~'^"~Mm ^ -'^i^* / + '^2(^ + |c|)' + fi2h\\v\\li . 

Note that (15. ip implies |c| < /i + He"^'^""'!'!!^!. Substituting this into (16. 8p yields 
(E3). ^ □ 

7. Energy estimate 
Recall that /C = 4c^ — 9^ — 6?7(-, a, c). Let 

Lemma 7.1 (energy estimate). Suppose that we are given bo G C^(M^) and 5 > 0. 
(Implicit constants below depend only on bo and 6.) Suppose v solves (14. ip and satisfies 
( im . and 5 < c < r/ien 

(7.1) < I - a + 4c^ - b{a)\\\e-^\^-Mli + ^Mk 

We remark that by integrating (17. ip over [0,T], 1 < T ^ and applying 

Corollary 12. 4^ we obtain 

(T.2) Wvh^^^^ni < ll^^o||H^ + ||d-4c^ + 6(a)||i|j|e-^l— 

Proof. We compute 

(9t£^(f ) = {ICv, dtv) - 3{v'^, dtv) + 4cc||f ||^2 - 3{{adar] + cdcr])v, v) 
= I + 11 + III + IV 
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Into I, we substitute fl4.ip . This gives 

I = (/Ct;, d,Xv) - Ac^{ICv, d^v) + {Kv, d^{bv)) - 3{ICv, d^) + (JCv, Fq) 
= lA + IB + IC + ID + IE 

We have lA = 0, while IB = —12c^{rix., v"^). For IC, numerous applications of integra- 
tion by parts gives 

IC = 2c^{b„ v^) + ^{h^, vl) - ^(6,,,, v^) - 3{7]b„ v^) + 3(r/^6, v^) , 
and hence 

lC = ?>h{a){r,.,,v^) + 0{h\\v\\l,). 

Note 

IE = {vXF\\) + {vXFi) . 

But since ICdaTj = 0, Kdci] = r], and {v,ri) = 0, we have {v,lCF\\) = 0. We estimate 
the second term to obtain 

|IE|</i||e-l--'^l^||Hi. 

Combining, we obtain 

I = (12c2 - 3b{a)){daV,v^) - 3{lCv,d.y) 
+ 0{h\\v\\]j, + h\\e^'\''-Mm) ■ 
Substituting (14. ip into II, we obtain: 

II = -3{v\ d^lCv) + I2c\v\ d^v) - 3{v\ d^ibv)) + d^) - 3{v\ Fq) 
In II, we keep only the first term and estimate the rest to obtain 

II = -3(^;^ d^lCv) + Oih\\v\\l, + ||e-^l^-'^l^;||^i HFoe^^l^"'^! . 

Note 

l|e'^'^-'^'Fo||^,i<|d-4c2-6(a)| + |c| + /i. 

Collecting, we obtain 

(7.3) \dt£\ <\-a + Ac'- b{a)\\\e~^\^-Mli + + l^DMk 

Note that in the addition of terms I and II, the terms ±(t;^, dxlCv) canceled, and in the 
addition of I and IV, the two 0(1) coefficients —3a and 12(? — 3b{a) were combined 
to give the smaller coefficient —3d + 12c^ — 36(a). 

Finally, we note that (15. ip implies |c| ^ h + ||e~''l'^~"u||^i. Substituting this into 
(17:31) yileds (EID. ^ □ 
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8. Proof of Theorem [2] 
It will be shown later that Theorem [2] follows from the following proposition. 

Proposition 8.1. Suppose we are given bo € C^(M^) and 6 > 0. (Implicit constants 
below depend only on b^ and 6). Suppose that we are further given G M, Cq > 0, 
K > 1, h > 0, and Vq satisfying f ll.6p . such that 

Let u{t) be the solution to (11. 2p with b{x, t) = b^lkx, ht) and initial data ri{-, a^, Co)+Vq. 
Then there exist a time T' > and trajectories a{t) and c{t) defined on [0,T'] such 
that a(0) = ao; c(0) = Cq and the following holds, with v u — ri{-, a, c): 

(1) On [0,T'], the orthogonality conditions (11. 6p hold. 

(2) Either c{T') = 5, c{T') = 5~\ or T' ~ h'^ . 
3 \\v\\l^ m < /^/^^/^ 

(4) \\e-'\''~''\v\\r2 ^1 < kH^/'^. 

(5) /o^' \d-Ac^ + b{a)\dt < k. 

(6) j^' \c-\cb'{a)\dt<K'^h. 

Proof. Recall our convention that implicit constants depend only on b^ and 5. By 
Lemma l3.ll and the continuity of the flow u{t) in H^., there exists some T" > on 
which a(t), c{t) can be defined so that (11.60 holds. Now take T" to be the maximal 
time on which a(t), c{t) can be defined so that (I1.6P holds. Let T' be first time 
< T' < T" such that c(T') = 5, c(T') = 5"^, T' = T", or uh~^ (whichever comes 
first). Here, < a; <^ 1 is a constant that will be chosen suitably small at the end of 
the proof (depending only upon implicit constants in the estimates, and hence only 
on bo and S). 

Remark 8.2. We will show that on [0,T'], we have 11^(^)11^1 < kH^^"^, and hence 
by Lemma 13.11 and the continuity of the u{t) flow, it must be the case that either 
c(T') = 5, c(T') = or ujh~^ (i.e. the case T' = T" does not arise). 

Let T, < T < T', be the maximal time such that 
(8.1) Ibllico ^1 < anh^/^ , 

where a is a suitably large constant related to the implicit constants in the estimates 
(and thus dependent only upon bo and 5 > 0). In fact a > 1 is taken to be 4 times 
the implicit constant in front of HfoUffi in the energy estimate (17.21) . 

Remark 8.3. We will show, assuming that (18. ip holds, that ||f HLjjyj/i'i ^ ^anh^^'^ and 
thus by continuity we must have T = T'. 
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In the remainder of the proof, we work on the time interval [0, T], and we are able 
to assume that the orthogonality conditions ( 11. 6p hold, 6 < c(t) < 6~^, and that ( 18. ip 
holds. We supress the a dependence in the estimates in (18 ■2p and (18. 3p below. 

By Lemma ISTT] (15. 2p . and (18. ip . (just using that ||e~''l^'~"k'||j|/i < Ht'll/fi) it follows 
that 

(8.2) \a - Ac + b{a)\ < Kh^^\ 

By (18. 2p . the hypothesis of the local virial estimate Lemma 16.31 is satisfied. Using 
flO) in (EJ]) (recall T = uh-^ < h'^), we obtain 

(8.3) \\e-'^''-MLf^^,,Hi < ^h'^' . 

Inserting (18. ip . (18. 2p . and (18. Sp into the energy estimate (17. 2p (recall T = uh~^), we 
obtain 

Provided h <q, k""^ and uj <^a 1, we obtain (recall ||t'o||//i ^ nh^^"^), we conclude that 
IK'Ulj^^j//! < completing the bootstrap, and demonstrating that T = T'. In 

particular, we have established items ([1]), (|2]), (|3]), (j4]) in the proposition statement. 
It remains to prove ([5]) and . By Lemma 15.11 (15. ip , 
rT 

/ |c-k6'(a)|rft</iTi/le-^l"-°li;L2 m + \\e-'\''-''\v\\% „i + T/i^ 
I 3 wi ~ II iii.[Q,j,jii^ II iiLjQ,j,jii^ 

(8.4) <K'h, 

establishing item (|6]). Similarly by Lemma [5.11 (15. 2p . we obtain item ([5]). □ 
The above proposition can be iterated to obtain: 

Corollary 8.4. Suppose we are given bo G C^(M^) and 6 > 0. (Implicit constants 
and the constant C below depend only on b^ and 6). Suppose that we are further given 
OoGM, Co > 0, /3 > 1, h > 0, and Vq satisfying (II. 6p . such that 

Let u(t) be the solution to (II. 2p with b{x, t) = bo{hx, ht) and initial data ri(-, gq, co)+'yo. 
Then there exist a time T' > and trajectories a{t) and c{t) defined on [0,T'] such 
that a(0) = flo; c(0) = Cq and the following holds, with v '= u — ri{-, a, c): 

(1) On [0,T'], the orthogonality conditions ( 11. 6p hold. 

(2) Either c{T') = 6, c{T') = 6~\ or T' ~ h'^ \ogh-^. 

(3) \\v\\,^^^,^nl<Ph^/'e^'\ 

(4) \\e-^\^-MLl^,,^Hi<Ph^'^e^^K 

(5) j'^' \a - 4c2 + 6(a) I dt < /3e^'^* 
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(6) J^' \c - lcb'{a)\ dt < /3^he^^'. 



Proof. Let K ^ Ihe the constant that appears in item ([3]) of Prop l8.ll and < w <^ 1 
be such that T' = ujh~^ in item of Prop. 18.11 Let kj = [5K^ for 1 < j < J, where 
J is such that K'^ ~ h~^^^. Let Ij denote the time interval Ij = [(j — l)ujh~^,jujh~^]. 
Apply Prop. 18. li on L with k = k,-. □ 



Now we complete the proof of Theorem |2l Recall that we are given bo G C^(M^) 
(5 > 0, ao G M, and Cq > 0. Let A{t), C(r), and be given as in Def. [H Let T', a(t) 
c{t) be as given in Cor. El Let d{t) = h'^A^ht) and c(t) = C{ht). Then 



a - 4c + b{a) 
c — |c6'(a) 











on < t < h-^T^. Then 

\a — d\(t) < [ \d — d\ds 

I (4c2 - 6(a)) - (4£2 - 6(5)) | (s) ds + / \d - 4c' + b{a) \ (s) ds 

< / \c — c\{s)ds + h / \a — d\ ds + /3'^e'' 
Jo Jo 



< 



'0 



By Gronwall's inequality, 
(8.5) 
Also, 



a-o|(t) < e^'^* ( / |c-c|(s)c/s + /32 ) . 



\c-m 



< 



In- 



< / \b'{a)-b'{d)\ds 



= Inc 


-ln£|(t) = 






/*!-- 


-b'{a) \ ds 


Jo c 




CM 





c c 
c c 



(s) ds 



< /i / \a — a\ ds + iS'^he 
Jo 

Combining, and applying Gronwall's inequality again, we obtain 

|c-£|(t)</3V^*. 
Substitution back into f l8.5p yields 

This completes the proof of Theorem [2l 
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Appendix A. Global well-posedness 

In this section, we prove that (11.21) is globally well-posed in H^. The local well- 
posedness (Prop. lA.ll below) is a consequence of the local smoothing and maximal 
function estimate of Kenig- Ponce- Vega [13] and the global well-posedness follows from 
the local well-posedness and the nearly conserved norm and Hamiltonian (Prop. 
lA. 21 below). A similar argument is given in Apx. A of [1] with an additional smallness 
assumption on b. This smallness assumption could be removed by scaling their result. 
However, for expository purposes we present a shorter proof here, which also imposes 
fewer hypotheses on b. 

In this section, we adopt the notation to mean L^qj^^ and CtH^ to mean 
C([0, T]; i7|), etc. The ordering of multiple norms is standard: for example, ||w||L2j;^|f = 

II lklU?f Wli- 

Proposition A.l (local well-posedness of (II. 2p in H^). Let X he the space of func- 
tions on [0, T] X M defined by the norm 

\\w\\x = IMl^l^ + |k||c,,[o.T]^^i 

Suppose that 

A = IIMIlIL-,, ^ + ||5xfe||L-;„,jLi < ^ 

and (j) G . Then there exists T = T{A, \\(f)\\{ji) < 1 and a solution u & X to (11.21) 
with initial data on [0,T]. This solution is the unique solution belonging to the 
function class X. Moreover, the data-to- solution map is Lipschitz continuous. 

Proof. Let U denote the linear flow (no potential) operator, a mapping from functions 
of X to functions of (x, t), defined by 

{U<P){x,t)=e-'''<P{x) = i-^e^-«'0(Orf^ 

Let / denote the Duhamel operator, a mapping from functions of (x, t) to functions 
of (x, t), defined by 

(J/)(x,t)= fe-('-''^^'f{.,t')dt'. 
Jo 

That is, iiw = U(f), then w solves the homogeneous initial- value problem dtw + d^w = 
with w{0,x) = 0(x). If w = If, then w solves the inhomogeneous initial-value 
problem dfiu + d^w = f with n(0, x) = 0. 

Kenig-Ponce-Vega [12l [13] establish the estimates 

(A.l) \m\crLi < Uhl 

(A.2) imui^ < mm 

(A.3) WdJfWcrLi < II/IIlil^ 

(A.4) \\dJf\\LlL^<\\f\\LlL^^ + \\dJhlL^^ 
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with implicit constants independent of < T < 1. In fact, f lA.ip is just the unitarity 
of U{t) on Ll, (1X21) is (2.12) in Cor. 2.9 in p], ( lAlal) is (3.7) in Theorem 3.5(ii) in 
|13j . and (IA.4|) is not exphcitly contained in [121 [IS]; but can be deduced from the 
above quoted estimates as follows. By the Christ-Kiselev lemma as stated and proved 
in Lemma 3 of Molinet-Ribaud |16| , it suffices to show that 



U{t-t')f{t') dt' 



< 



T 2 T oo 



By first applying flA.2p and then the dual to the local smoothing estimate \\dxU 011 00^2 



< 



11011^2 (Lemma 2.1 in [12]), we obtain 



< 



U{t-t')f{t')dt' 
d. 



T 



t')f{t')dt' 


+ 


[ 






Jo 



U{-t')dJ{t') dt' 



Ll 



as claimed. 

Let $ be the mapping 

(A.5) 



$(«;) = U(t) + dj{w^ - hw) 



We seek a fixed point $(m) = m in some ball in the space X. To control inhimo- 
geneities, we need the following four estimates, which are consequences of Holder's 
inequality: 



mhu)\\^.^.^<Ty\\\dML^Ll\\u\\Ll 
\M\lILI < T^''^\\'b\\LlL^\\u\\L^Ll 

I 2|| ^ rrl/2\\ II II II 



(A.6) 
(A.7) 
(A.8) 
(A.9) 

We prove flA.6l) . 

\\dx{bu)\\LiLl < \\dxb\\L2Ll\\u\\LlL^ + \\b\\LlL^\\dxU\\LUl 

< \\dxb\\L2,Ll\\u\\LlL^ + \\b\\LlL^\\dxU\\L2,Ll 

< T^^'^\\dxb\\L^Ll\\u\\LlLS? +T^^'^\\b\\L'iL^\\dxU\\L^Ll 



which is flA.6l) . The other estimates, flA.7p . flA.Sp . flA.91) are proved similarly. 
By (ED, dED, 



(A.IO) 



\^(w) 



|l2L- < + ||(W -bw)\\LiLl + \\dx{w -MILiL2 
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By dH]), (jO), 



(A.ll) mw)\\L^Ll < + U^' - MlkiL^ 

Applying dx to flA.SP and estimating with flA.ip . flA.3l) . 



(A.12) \\dMy^)\\L¥^ < \\d.<P\\Li + \mw' - MILiL^ 

Combining (lA.lOl) . (lA.llI) . (1A.12I) . and bounding the right-hand sides using (lA.Gp 
(OTTj) . flXSj) . fOTOD . we obtain 

(A.13) ||*H||x < + CT'/\A\\w\\x + \\w\ 



Xj 



Let B = 2C||0||hi, and consider Xb = {w eX \ \\w\\x < B} and T < ^C'^ mm{A-^, B~ 
Then (lAl3l) imphes that $ : ^ X^. 

We similarly establish that $ is a contraction on Xb, which completes the proof. □ 

Proposition A. 2 (global well-posedness of fll.2p in H^). Suppose that b E C^(M}^^) 
and satisfies the following. Suppose that for every unit-sized time interval I , we have 

(the bound need not be uniform with respect to all time intervals). Also suppose that 
for all t, 

|IW)IU- < oo, IIW)IU- < oo- 
Let (j) G . Then the local solution to fll.2p with initial data given by Prop. 
\A.1\ extends to a global solution with 



where 7(5) is given by 



lif) = \\bx{s)\\L^^^^L^ds. 

Proof. Let P{t) = \\u{t)\\'j^2 (the momentum) and recall the definition fILSp of H, the 
Hamiltonian. Direct computation shows that 

dtP = J bxU^ dx , dtH ^ ^ J ^tU^ dx . 

Then |P'(t)| < i{t)P{t), and hence dt[e-^^^^ P{t)] < 0. From this, we conclude that 

P{t) < e^Wp(O) . 

In addition, we have 

\H\t)\ < \\hmL^P{t) < ||fct(t)||Lgoe^Wp(0) 
Hence ^ 

if(t) < ff(0) + P(0) ! \\bt{s)\\L^e^'^'Us 
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By the Gagliardo-Nirenberg inequality ||m|||3 < and the Peter- Paul 

inequality, we have ||m|||3 < 111^x11^2 + Cll^ll^^''^- Hence 

hx.|li|<C||n||if + ||6(t)|Ug=||n(t)||i.+i7(t) 
When combined with the inequalities for H{t) and P{t), this gives the conclusion. □ 
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